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Abstract 

A meson can decay to K^ir through doubly Cabibbo-suppressed decay or via 
mixing. With 46.2 fb~^ of integrated luminosity collected by Belle, we have measured the time 
integrated rate of the wrong-sign process D'^ — > K+tt" relative to that of the Cabibbo-favored 
process K-'k+ to be Rws = (o.372 ± 0.0251^;^?^) % (preliminary). The D^-D° mixing 

parameters can be derived from the time distribution of the wrong-sign process. 

PACS numbers: 13.25.Ft, 12.15.Ff, 14.40.Lb 
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I. INTRODUCTION 



The Standard Model predicts that D^—D^ mixing is small, with many estimates for the 
mixing parameters x = AM/T and y = Ar/2r at the ^ 10~^ level, although long-distance 
effects in the Standard Model may raise both parameters to ~ 10^^ ^. New physics 
effects may enhance x, but are not expected to affect y. Consequently, a discovery of large 
X with small y would imply new physics. 

Experimentally, mixing can be identified by the study of the wrong-sign (WS) process 
D° K^n^ , which can proceed directly through doubly Cabibbo-suppressed decay (DCSD) 
or through mixing followed by Cabibbo-favored decay (CFD). The decay of an initial D*^ — >■ 
D°7r+ yields a charged "slow" pion, denoted 7r+, whose sign can be used to tag the initial 
D as either D° or D^. 

The differential WS rate relative to the right-sign process is |^ 

T{D\t) K+vr-) 

rws{t) ^ 



r(L)0 K-Ti^ 



Rd + 4^y'Tt + h^x'^ + y'^)VH'' 



(1) 



where Rd is the relative rate of DCSD to CFD. We use the convention of CLEO [0, with 
y' = y cos (5 — X sin (5 and x' = x cos 5 + y sin S, where S is the relative strong phase between 
the DCSD and CFD channels. A fit to the rws{t) distribution in Belle data, to obtain the 
mixing parameters, is underway. In this paper, we present a preliminary measurement of 
the time-integrated wrong-sign rate, 

Rws = Rd + \fR^y' + \ {x^ + y") , (2) 

based on the same event reconstruction and background-handling techniques used for the 
full analysis. 



II. EXPERIMENTAL APPARATUS AND DATA SAMPLE 

The analysis is based on data accumulated by the Belle detector at KEKB [Q, an asym- 
metric collider with 8 GeV electron and 3.5 GeV positron storage rings. The data set 
corresponds to an integrated luminosity of 46.2 fb~^. 

Belle 10] is a general-purpose detector based on a 1.5 T superconducting solenoid. Charged 
tracks are reconstructed using a 50-layer Central Drift Chamber (CDC) and a 3-layer double 
sided Silicon Vertex Detector (SVD). Particle identification for charged kaons and pions is 
performed by combining information from the CDC {dE/dx), a set of time-of- flight counters 
(TOF) and a set of aerogel Cerenkov counters (ACC). The combined result of the three 
systems provides K/tt separation up to 3.5 GeV. 



III. CANDIDATE RECONSTRUCTION 



For D'^ K^Tv decay, we use oppositely charged tracks where both tracks have good 
SVD hits. Charged K and vr mesons are required to be positively identified, with efficiencies 
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88.0% and 88.5% respectively, and fake rates of 8.5% (vr fakes K) and 8.8% {K fakes vr). 
D*^ candidates are then reconstructed by combining the candidate with another charged 
track, 7r+. The reconstructed D*~^ momentum in the center of mass frame P*{D*) is required 
to be greater than 2.5 GeV to ehminate BB events. Vertex fits are performed, requiring 
the two tracks forming the D'^ candidate to originate from a common vertex, the flight 
direction to be consistent with a particle originating from the interaction point, and the tt^ 
to originate from the Z)° production vertex. We accept candidates with good vertex quality. 
The value of mass M and the energy released in D*~^ decay Q = M{K~^ ,tt~ ,71^) — 
M{K+,TT-) - Af(7r+) are then obtained. 

The right-sign decay D*^ D^^^t •> ~^ K~tt^ is reconstructed in the same way. Here 
and throughout this paper, the inclusion of charge conjugate modes is implied. 

IV. BACKGROUND DETERMINATION 

The information in the signal region alone is insufficient to determine the signal time 
distribution, because the background level is comparable to the signal in the signal box, and 
the various backgrounds have different time distributions. We therefore divide the back- 
ground into categories, and determine the population of each using M, Q information. To 
obtain the signal yield as well as the level of each background, we perform a two-dimensional 
(2D) fit to (M, Q). By studying a large number of continuum Monte-Carlo (MC) events, we 
categorize the backgrounds as follows: 

• —>■ K^Ti^ with double misidentification {K~ 7r^,7r^ —>■ K^). 

• —>■ K^K~ with single misidentification {K~ ^ tt~). 

• TT+TT^ with single misidentification (tt^ K^). 

• K^TT^ combines with a random slow 7r+ to form a fake D*~^. 

• > 3 body decay, with two daughters identified as K, vr. 

• D+, decay, with two daughters identified as K, n. 

• Pure combinatorial background. 

The — > K~7i~^ double misidentification background can be rejected by kinematic 
requirements. We evaluate the D° mass under swapped mass assignment Mf^p = Mk^-k,-k^k- 
If Mflip falls within 28 MeV (~ 4cr) of the nominal mass, the K^tt^ candidate 

is rejected. For the K^K~ ,tt^tt~ single misidentification background, the change 

from the hypothesis that a charged track is a 7i{K) to a K{n) in the 2 body decay 
daughters results in an increase(decrease) in the reconstructed mass M of at least 60 
MeV. This behavior has been studied in both data and MC. So we choose the mass 
window 1.81 ~ 1.91 GeV for histogram making and fitting to exclude these two backgrounds. 

Each of the remaining four types of background is parametrized and fitted with a 2D 
function to obtain its shape, as shown in Figs |l]-^ The backgrounds in the right-sign 
sample are categorized and fitted in a similar way. They are grouped into three types: 

• —>■ K^TT^ combines with a random slow tt^ to form a fake D*~^. 
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FIG. 1: The Q distribution for MC events (points) and the fit function (histogram) for the random 
slow tt"*" background. The M shape is fixed to that of the right-sign signal in the data. 




3hody Mass GeV 3hody Q MeV 



FIG. 2: M (left) and Q (right) distributions for MC events (points) and the fit functions (his- 
tograms) for the > 3 body background. 
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FIG. 3: M (left) and Q (right) distributions for MC events (points) and the fit functions (his- 
tograms) for the Df,D^ background. 




FIG. 4: M (left) and Q (right) distributions for MC events (points) and the fit functions (his- 
tograms) for the combinatorial background. 
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• Charmed ground state meson > 3 body decay, with two daughters identified as K, vr 
("D >3 body" background). 



• Pure combinatorial background. 



V. RESULTS 

We perform a 2D fit to the (M, Q) distribution for the right-sign data, with the nor- 
mafization of each background aUowed to fioat in the fit. The signal is represented by a 
double Gaussian in M and a bifurcated Student's t function in Q. The projections onto 
M and Q are shown in Fig. ^ To perform the 2D fit for the wrong-sign data, we fix the 
wrong-sign signal shape to the shape of the right-sign signal, and fioat the normalization 
of the signal and each type of background, but fix the relative normalization between the 
D^, and > 3 body backgrounds. The projections of the wrong-sign data and the fit 
results are shown in Fig. ^. The data are the circles with error bars, and there are 20 bins in 
M (1.81 ~ 1.91 GeV) and 160 bins in Q (0 ~ 20 MeV). We find 120795 ± 371 ^ K-7i+ 
and 450 ±31 K+tt- events. The ratio Rws is then calculated to be (0.372 ± 0.025) %. 




1.825 1.85 1.875 1.9 5 10 15 20 

Mass (GeV) Q (MeV) 



FIG. 5: Projections of M (left) and Q (right) for the right-sign data (points) and the fit functions 
(histograms), for the region 1.81 < M < 1.91 GeV and < Q < 20 MeV. Note the logarithmic 
scale. 

Many systematics cancel because Rws is a ratio between two decay modes with similar 
kinematics. The dominant systematic errors stem from potential imperfect modeling of the 
shapes for our backgrounds. We vary selection criteria {K and vr identification cuts, vertexing 
criteria, the P*{D*) cut) and use different background parameterizations, and then repeat 
the fit on each category of MC background as well as on the right-sign and wrong-sign data. 
The variation in the result is taken as a measure of the systematic error. We also vary the 
background shape parameters according to their errors and calculate the resulting error on 



8 




Mass (GeV) Q (MeV) 



FIG. 6: Projections of M (left) and Q (right) for the wrong-sign data (points) and the fit 
functions (histograms), within a 3a window in the complementary variable (5.27 < Q < 6.47 MeV 
and 1.8445 < M < 1.8845 GeV respectively). The signal contribution is shaded yellow. 



TABLE I: Summary of systematic errors on Rws 



Source Systematic error (%) 

Kaon identification io 0023 

Pion identification -aoo48 

Y . n. +o:oo42 

veiT^ex nt -0.0051 

P*{D*) cutoff t'o fiTr 

Background shapes ^o 0026 

Background parametrizations ±0.0011 

T^T—j +0:009 
lotai -0.014 



Rws- The correlations among these background shape parameters are considered either in 
the variation or in the Rws error calculation. The results are summarized in Table |. All 
systematic errors are combined in quadrature. 

In summary, we have measured the ratio Rws of the rate of the wrong-sign process 
K^7T~ relative to the right-sign process K^tt'^ to be 

Rws = (0.372 ± 0.025(stat);°:°°^(syst)) %. (3) 

This result is preliminary. It is compared with other experimental results |^, ^ ^ in Fig |^. 
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FIG. 7: Various experimental results for Rws- 
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